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Abstract 

We generalize the model of solid inflation to an anisotropic cosmic solid. Barring 
fine tunings, the observed isotropy of the cosmological background and of the scalar 
two-point function isolate the icosahedral group as the only possible symmetry group of 
such a solid. In such a case, higher-point correlation functions—starting with the three- 
point one—are naturally maximally anisotropic, which makes the standard detection 
strategies highly inefficient and calls for a dedicated analysis of CMB data. The tensor 
two-point function can also be highly anisotropic, but only in the presence of sizable 
higher-derivative couplings. 


1 Introduction 

The universe appears isotropic on large scales, and it is thus natural to assume that whatever 
it was that fueled primordial inflation, it was an isotropic system. It is interesting, however, 
to analyze critically this assumption. Observations tell us that the cosmological background 
and the spectrum of scalar perturbations are isotropic, but they do not tell us anything about 
higher-point correlation functions or about tensor modes, for the simple reason that we have 
not detected them yet. 

This raises the following question: can one have a physical system driving inflation whose 
dynamics are intrinsically anisotropic—perhaps maximally so—but that nevertheless yields 
an isotropic background and an isotropic scalar spectrum of perturbations? We are not 
interested in systems for which one can tune coefficients in the Lagrangian in order to achieve 
the desired degree of isotropy, but rather in systems whose symmetries are so powerful 
as to enforce such an isotropy, leaving open the possibility of anisotropic signals for other 
observables. 

To make the discussion more concrete, let’s consider the cubic group. This is the discrete 
subgroup of rotations that maps a cubic lattice into itself. Calling x, y, z the lattice’s 
preferred directions, the cubic group is simply the set of permutations of x, y, z as well as 
single-axis inversions x —)■ —x, etc. Barring fine tunings, the dynamics of an homogenous 
system with this symmetry group—such as a cubic crystal in the continuum limit—in general 
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will not be isotropic. However, certain observables are forced to be. In particular, because of 
cubic symmetry, any two-index tensor associated with the lattice must take the form 

oc + y^'y^ + = 6^^ . (1.1) 


On the other hand, with more indices there are structures that are invariant under the 
cubic group only, and can in principle lead to observable anisotropies. For instance, at the 
four-index level, the tensor structure 


X x^x X + y y^y y + z zz z z 


( 1 . 2 ) 


is invariant under the cubic group, but cannot be rewritten in terms of Kronecker deltas only. 

Solids are natural candidates for considering discrete subgroup of rotations, and for this 
reason we will elaborate on the above ideas in the context of solid inflation [ 1 ] , which we 
review in the next section. In solid inflation, inflation is driven by a solid’s stress-energy 
tensor . For the background evolution to be isotropic, one needs an isotropic T*-^ on the 
ground state of the solid. However, in general is invariant only under the symmetries 
of the solid under consideration, which restricts the number of possible symmetry groups 
to those whose invariant two-index tensors are accidentally isotropic. As we saw above, the 
cubic group has this property. 

Moving away from the background, the constraints become more severe. The fluctuating 
degrees of freedom are the solid’s phonons, which can be parametrized by a 3-vector field ^{x), 
and the metric perturbations. To discuss possible anisotropies of scalar correlation functions, 
it is sufficient to focus on the phonons; the longitudinal one mixes with the scalar metric 
perturbations, and so any anisotropies in its dynamics will be reflected in scalar correlation 
functions. In particular, the two-point function is determined by the phonons’ quadratic 
Lagrangian, which takes the general form 

C 2 = AijTT^^ + , (1.3) 

where and tensors that are invariant under the symmetry group of the solid. 

We see that for the scalar two-point function to be isotropic, one also needs the invariant 
four-index tensors to be isotropic. As we saw above, the cubic group does not pass this test. 

We will show in sect. 3 that the only discrete subgroup of rotations with the above 
properties is the icosahedral group—the symmetry group of an icosahedron. The natural 
question now is which observables are going to exhibit the anisotropies associated with such 
an icosahedron: there should be many preferred directions in the sky (20, 30, or 12, depending 
on whether one counts the faces, edges, or vertices), which should show up in correlation 
functions. The question can be approached once again in terms of invariant tensors. The 
scalar three-point function is determined by the cubic Lagrangian for the phonons, which 
now can involve a six-index tensor: 


C 3 D • (1.4) 

We will show that the icosahedral group allows for the anisotropic invariant tensor 

OC 2(7 + 2) + (7 + 1) + ...) + (^^ijkl^mn^mi -1 ^ ^ 
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where 7 is the golden ratio, the six-index and four-index deltas are nonzero only if all their 
indices take the same value, and i + 1 and i — 1 are to be interpreted modulo 3, that is, 
3 -|- 1 —)■ 1 and 1 — 1 —)■ 3. This makes the three-point function potentially highly anisotropic. 
In fact, we will show that there is a choice of Lagrangian coefficients for which the three-point 
function is completely anisotropic, in the sense that it has exactly vanishing overlap with any 
three-point function template associated with isotropic models. 

Anisotropies can show up in the tensor spectrum as well. The reason is that the quadratic 
Lagrangian for tensor modes also involves a six-index invariant tensor, 

^ daiudamn. (1.5) 

However, we will see that such six-index tensor can receive anisotropic contributions only 
from higher-derivative terms in the Lagrangian. It is consistent within our effective theory 
to assume that these are so large as to yield order-one anisotropies in the tensor spectrum, 
and in the examples we considered we found no indication that this disrupts the technical 
naturalness of the effective theory, but we have not investigated the question systematically. 
We leave this for future work. Notice that a strong anisotropy in the tensor spectrum can in 
principle reconcile tensions between a large tensor signal in a small patch of the sky and little 
or no signal in a whole-sky average, like the original—then evaporated—BICEP2/Planck 
tension. 

Finally, we should emphasize that when we talk about ‘solids’ we do not mean systems 
with an underlying crystal structure, but rather continuous, homogeneous solid media, which 
can be more symmetric than crystals in the continuum limit. For instance, there is no crystal 
with icosahedral symmetry group, but it is perfectly consistent to assume such a symmetry 
for a continuous medium (in fact, there are ^tiosz-crystals with icosahedral symmetry [ 2 ].) 
Even though the solids of everyday life are not homogeneous at microscopic scales, there is 
no a priori reason why there couldn’t exist (perhaps strongly coupled) held theories that at 
hnite density exhibit perfectly homogeneous solid-like states. If one is uncomfortable with 
such an assumption, one can regard our inhationary model simply as a system of three scalar 
helds with certain symmetries. As we now review, the low-energy ehective held theory is the 
same, which makes the diherence between the two viewpoints unsubstantial. 

2 Solid inflation 

From an ehective held theory standpoint, the mechanical deformations of an homogeneous 
solid can be described in terms of three scalar helds 0^(x) (/ = 1, 2, 3) [3], whose expectation 
values in the ground state of the solid are 

( 2 . 1 ) 

and whose Lagrangian is invariant under the shift symmetries 

(j)^ ^ (j)^ + a\ = const ( 2 . 2 ) 

(see also [4, 5] for alternative approaches.) The 0'^’s can be regarded as the comoving 
coordinates of the solid’s volume elements. By PoincarA and shift-invariance, to lowest order 
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in derivatives the Lagrangian must take the form 

C = , (2.3) 

where F is an a-priori generic function, determined by the solid’s equation of state. For a 
solid with symmetry group G C 50(3), one also demands that the Lagrangian be invariant 
under the internal rotations 


O'jeG. 


( 2 . 4 ) 


which restricts the form of F. 

For instance, in the case in which G is the full 50(3)—which is the case extensively 
studied in [1 ]— F can only depend on three invariants, e.g. 

[B] , [B^] , [B^] , (2.5) 

where the square brackets denote the trace of the matrix within. Any other rotationally 
invariant function of B^'^ can be expressed in terms of these, e.g. 

detB = l{[Bf-3[B][B^] + 2[B^]) . (2.6) 

Upon minimally coupling the solid to gravity, the form of F is restricted further by 
demanding that the solid be able to drive near exponential inflation. In order for that 
to happen, one needs a solid that can be stretched by a large exponential factor without 
changing too much its physical properties, such as its energy density. Such a behavior is of 
course unlike that of any standard solid we know of, but it can be achieved by imposing an 
approximate internal scale invariance [1] 

(2.7) 


Focusing again on the 50(3) invariant case, to implement this symmetry it is useful to 
organize the three invariants (2.5) as 


X = [S] , 




[Bf 


( 2 . 8 ) 


X depends on the overall normalization of i?, but Y and Z do not, and as a result Y and Z 
are invariant under the transformation (2.7). The requirement of approximate scale invariance 
thus translates into a weak dependence of F{X, Y, Z) on X. In particular, by evaluating the 
solid’s stress-energy tensor on an FRW background, one hnds [1] 


p=-F, p = F-^Fx {X = 3/a^), (2.9) 


which yields the slow-roll parameter 


H _ 3 Fx _ dlogF 
IP ~ ~ dlogX ■ 


( 2 . 10 ) 


Of course “slow-roll” here is a bad characterization, because nothing is rolling, slowly or 
otherwise: the background conhgurations for our 0^’s only depend on the spatial coordinates. 
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But slow-roll parameters like e and the higher order ones can also be dehned geometrically, 
without any reference to rolling helds, purely in terms of the time-dependence of the Hubble 
scale H. We will adopt these geometric dehnitions—as we did above for e—and still use the 
standard slow-roll nomenclature. 

In the presence of perturbations, the scalars and the metric are 

(j)^ = x^ + 7 t\x) , (t) + 6g^^{x) , (2.11) 

where = diag(—1, a^, a^, a^) is the standard FRW metric. At distances much shorter 
than the Hubble radius, one can neglect the metric perturbations and identify 7f(x) with the 
phonon held. Expanding the solid Lagrangian to quadratic order, one hnds two parameters 
Cl and ct characterizing the longitudinal and transverse phonon propagation speed. These 
are determined by certain derivatives of F, evaluated on the background conhguration. In 
particular, one hnds the universal exact relation [1] 

4 = 1(1 + 4 - |e +, (2.12) 

where g = ejeH is the second slow-roll parameter. 

To study cosmological perturbations and compute their correlation functions, it is conve¬ 
nient to decompose the metric in an ADM fashion, 

ds^ = —N‘^dt^ + hij{dx^ + N^dt) [dx^ + N^dt) , (2-13) 

and choose spatially hat slice gauge (SFSG), 

hij = a(t)^exp (7y) , d^ij = 7** = 0 . (2.14) 


The curvature perturbation ^ can be dehned in a gauge-invariant fashion, and in the above 
gauge it is related to the phonon held vf by 

C = |V-7f. (2.15) 

Then, using standard cosmological perturbation theory, one can compute the correlation 
functions for scalar and tensor modes. At the two-point function level, the relevant observables 
are the scalar tilt, the tensor tilt, and the tensor-to-scalar ratio: 


ns — 1 ~ 2 e4 — 5s — 7 (2.16) 

nr - 1 ^ 2 c\e (2.17) 

r~16e4, (2.18) 

where s monitors the time-dependence of c^, s = clIclH. Particularly unusual predictions 
are the positivity of the tensor tilt—which would usually require a violation of the null energy 
condition—and the strong suppression of the tensor-to-scalar ratio in the slow sound speed 
limit, a factor of c\ stronger than for standard single-held cases. 

Expanding further the solid Lagrangian to cubic order, one hnds that at leading order in 
slow-roll the phonon self-interactions take the form 


TF) = ]Vfpa(t)^i7^^| A(5.7r*)3_ _ ^diTT^djTr’^dj'ir’^+ ^dj7r^dj7r'^dkTT^ 


(2.19) 
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The gravitational corrections to this are suppressed both in the de-mixing regime, k 3> 
and in the strong mixing one, k and one can argue that the cubic Lagrangian 

above is all one needs to compute the three-point function of curvature perturbations [1]. 
To leading order is slow roll, the result is 


(C(fcl)C(fc2)C(fc3)> ^ {2T:f5\k, + k 2 + h) 


3 Fy 1 
32 F 


X 


Q{ki, k2, k3)U{ki, k2, ks) 

U3U3U3 


where 


Q{ki,k2,k3) = ^kik2k3 - ^ 


81 


27 




+ k 


2/^3 


{h ■ hf 


+ k 


(fci • k 2 
kik 2 


2 {ki ■ k2){k2 ■ k3){k3 ■ ki] 

3 


kik2k3 


and 


U{ki,k 2 ,k 3 ) = 


3{k^^ + kl + kl) +20klklk. 


kik 2 k 3 {ki + k 2 + 

-|-18(fc^A;2^3 + -|- kik2kf) + 12{k\k2 + k^k^ -|- klkf) 

9 {k\k 2 + Sperms) -|- 12{k{k2 + Sperms) 

-|- 18 (fc^/c 2 ^ 3 -I-Sperms) I . 


Assuming Fy ~ F, this has a potentially huge /nl, 

1941SFy 1 


/nl — 


13122 F eel 


ec 


2 ’ 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


but its most peculiar feature is probably its ‘shape’ [6]— in particular, its purely quadrupolar 
angular dependence in the squeezed limit ks < kix- 


(CCC) « 

where 6 is the angle between ki and /C3. 


(1 — 3cos^ 9) 
TTp 


(2.24) 


3 Hunting for the right symmetry group 

We now want to generalize all of the above to a more general solid, invariant only under a 
discrete subgroup of rotations, which nonetheless features the desiderata identified in the 
Introduction: an isotropic background stress-tensor, and an isotropic quadratic Lagrangian 
for the phonons. As we saw, at the mathematical level these requirements are equivalent 
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to demanding that all invariant two-index and four-index tensors be fully isotropic for the 
symmetry group in question. 

There is an infinite number of discrete subgroups of 50(3), divided into two main classes: 
the crystallographic point groups and non-crystallographic ones. Let’s start with the former 
class. A crystallographic point group is the symmetry group of a crystal system that can hll 
all of space. This means that the group has to map all the lattice points into one another, 
which is a stronger requirement than being simply a subgroup of rotations. Since there is only 
a finite number of crystal systems—triclinic, monoclinic, orthorhombic, tetragonal, trigonal, 
hexagonal, and cubic—there is only a hnite number of crystallographic point groups. Except 
for the hexagonal one, all crystal systems can be dehned in terms of their three primitive 
lattice basis vector, let’s call them a, b, and c. Then, the two-index tensor 

+ UV + cV (3.1) 

is invariant under the corresponding crystallographic point group. However, this tensor is not 
invariant under general 50(3) rotations unless a, b, and c all have the same length and are 
all orthogonal to one another. So, only the cubic crystal survives. Still, as already pointed 
out in the Introduction, the cubic group fails our test at the four-index level, because the 
tensor structure 

a^a^a}a^ - 1 - bVb%'^ + cVc'c"* (3.2) 

is invariant under the cubic group but not under general S0{3) rotations. We thus reach the 
conclusion that no crystallographic point group can meet both of our criteria 

The non-crystallographic point groups are the icosahedral group, the inhnitely many Cn 
groups (u-fold rotations about a given axis), and the extensions of Cn that include some kind 
of reflection. In the two last cases, already at the two-index level we can easily construct 
invariant tensors that are not 50(3) invariant: for instance, the projector onto the plane 
perpendicular to the rotation axis. 

So, all our bets are on the icosahedral group—the symmetry group of the icosahedron. 
The icosahedron has 20 triangular faces, 30 edges, and 12 vertices, and there are 60 proper 
rotations that maps it into itself. Following [7], we orient our cartesian axes so that the 
icosahedron has two sides parallel to x, two parallel to y, and two parallel to z, as depicted 
in fig. 1. In this case, the coordinates of the vertices are (up to an overall rescaling) 

(± 7 , ±1,0) , (0, ± 7 , ±1) , (±1,0, ± 7 ) , (3.4) 

where 7 = is the golden ratio. 

We can hnd the invariant tensors in the following way. By dehnition, the two-index 
invariant tensors should satisfy 

= I\ P (3.5) 

^The only possible exception to this argument is the hexagonal crystal, some of whose links are not 
primitive lattice vectors but rather suitable linear combinations thereof. Still, one can easily show that certain 
two-index tensors that are invariant under the hexagonal group are not SO{3) invariant, e.g. 

^b\b{ + b^bi + bibi, (3.3) 

where 61 = (1,0,0), 62 = (1/2, ^3/2, 0), 63 = (-1/2, V3/2,0). 
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Figure 1: The relative orientation between our coordinate system and the icosahedron discussed in 
the text. 


for each rotation matrix I that belongs to the icosahedral gronp (we refer the reader to 
[7] for the explicit form of the rotation matrices). Given that there are 60 elements in the 
icosahedral gronp and 9 = 3x3 entries in T*b eq. (3.5) can be interpreted as a system of 
60 X 9 linear equations for the entries of T*b The solutions are ^ 

Til ^ y22 ^ T^j ^ Q 

That is, all two-index tensors that are invariant nnder the icosahedral group have the form 

oc , (3.7) 


and are therefore fnlly isotropic. 

We can apply the same logic to the fonr-index invariant tensors: 


rjiijlm _ _ ji jj jl jm rpabcd 


(3.8) 


In this case, we have 60 x 81 linear equations for the 81 = 3^ entries of and the solutions 

are 


211122 _ ^11133 _ 212211 _ 2^2233 _ 2^3311 _ 2^3322 

21I212 _ 2^1313 _ 2^2121 _ 2^2323 _ 2^3131 _ 2"3232 

2^1221 _ 2^1331 _ 2^2112 _ 2^2332 _ 2^3113 _ 2^3223 

21IIII _ 2^2222 _ 2^3333 _ 2^1122 2^1212 2^1221 

all Other = 0 . (3.9) 

These conditions can be rewritten compactly using only Kronecker deltas, which shows that 
all four-index invariant tensors are fnlly isotropic as well: 

Tijlm _ ^ _|_ (J 

^We used Mathematica to solve this linear system and those that follow. 
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(3.10) 









for arbitrary A, B, and C. 

Finally, let us consider the six-index invariant tensors. Following the same logic, we now 
have 60 X 3® = 60 X 729 equations 


rpijklmn _ rpfijklmn _ ji jj ^ jk jl^ jm j 


(3.11) 


On top of the isotropic solutions, schematically of the form T ~ 666, we hnd an anisotropic 
one: 


= 2(7 + 2 ) ^ -^ -^ 

where the dots stand for all other combinations of four and two indices out of six, the delta 
tensors with more than two indices are 1 only if all those indices take the same value, and 
i + 1 and i — 1 are to be interpreted modulo 3, that is 3 -|- 1 = 1 and 1 — 1 = 3. 

It is worth mentioning that one could have derived the invariant tensors above in a perhaps 
more intuitive fashion, by using as building blocks the 12 vectors Fa (a = 1 ,..., 12) that 
dehne the icosahedron’s vertices. Clearly, by taking suitable tensor products and summing 
over all the vertices, e.g. 

= (3.13) 

a 

one gets tensors that are invariant under the icosahedral group. However, it is not obvious 
that one can get all the invariant tensors in this way. Our brute-force analysis above settles 
the question (and the answer is ‘yes’, at least up to the six-index level, if one includes tensor 
products of lower-order tensors as well.) 

In conclusion, the icosahedral group has exactly the properties that we are after; all its 
two-index and four-index invariant tensors are isotropic, whereas its six-index ones are not. 
And it is the only subgroup of S'0(3) with these properties. From now on we will thus focus 
on a variant of solid inflation with icosahedral symmetry, which we dub ‘icosahedral inflation’. 
As already emphasized, the anisotropy of the six-index invariant tensors translates into an 
anisotropy of the scalar three-point function. Our goal now is to compute such a three-point 
function. 


4 Cubic Lagrangian of scalar modes 

An anisotropic invariant six-index tensor can induce anisotropies in the scalar three-point 
function through the trilinear phonon interaction 

£3 oc Tll^^^'^diTijdkTridmTrn ■ (4.1) 

However, to figure out the most general structure of the cubic Lagrangian compatible with 
our symmetries requires some work. In the S'O(3)-invariant version of solid inflation, this 
task was straightforward; the full solid’s Lagrangian only depends on the three invariants 
(2.5), each of which can be expanded in perturbations about the background solution, up to 
any desired order. In icosahedral inflation, we face the problem of classifying the allowed 
invariants of . Since B^'^ starts at zeroth order in if, 

B 1 + dn + dirdir , (4.2) 
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to expand the Lagrangian to any given order in ^ —cnbic, in onr case—we need to consider 
all orders in . However, at high orders, in principle we have to include more and more 
invariants, 

, (4.3) 

where T is a generic 2n-index tensor with icosahedral symmetry. We are not aware of any 
simplifying property of the icosahedral group analogous to the SO (3) statement that an 
arbitrary invariant of B^^ can be written as a non-linear function of the three fundamental 
invariants (2.5). Clearly, the number of independent invariants cannot be more than the 
number of independent components of B^^ —six—but using the individual components of 
B^'^ would make our computations messy and unreadable. 

To get around this problem, we can work directly with the fluctuation of B^^ about its 
background, but, as we will see below, we will have to be careful about the non-linearly 
realized symmetries In SFSG gauge, our building block up to cubic order is 

-N’^dkclS)+h'^^dkcl>^d^clS 
1 . . 1 

~ ^ -I—^ (cIV-^ -I- 5%*) H— -dkTT^d^TT^ — (tt* — N"") (tt^ — W) 

+ {r - W) N^dkTT^ + {tt^ - W) N’^dkTT^ + 2SN {r - W) (tt^' - W) -f- • ■ ■ 

= (4-4) 


where we dehned tt®-^ as the fluctuating part of B^\ and we stopped differentiating between 
the internal indices and the spacial i,j, ■ ■ ■ ones (the reason is that the background 

(0^) = breaks spatial rotations and internal ones down to the diagonal combination.) 
Were we to expand the full non-linear action up to cubic order, we would have 


F(H®^') = Fo + 


dF 


dB^i 


71 




d^F 


0 a2 ^ 2! dB^iB^^ 


0 


d^F 


3! dB^W^^B^ 




+ ■■■ , (4.5) 


where the subscript zeros mean ‘evaluated on the background’. By the background Friedmann 
equations, the hrst derivative of F with respect to B^^ can be related to e [1]: 


dF 

dB^i 0 


e p^Fo (5ij . 


( 4 . 6 ) 


This result was derived for the original solid inflation model assuming SO{3) invariance, but 
in the Appendix we prove that it holds for our icosahedral inflation case as well. The higher 
derivatives of F do not enter the Friedmann equations, and therefore cannot be related simply 
to other background quantities. As we will see below, they do obey constraints coming from 
the non-linearly realized symmetries, but for the moment we can just parametrize them as 
the most general icosahedral-invariant tensors with the right index-permutation symmetries 

^In the standard effective field theory of inflation [8], it is straightforward to write the action directly in 
terms of the metric perturbations in unitary gauge, by using for instance + 1 . As emphasized in 

[1], in solid inflation the analogous variable in unitary gauge would be Sg"^^ = g^^ — but this, unlike the 

full g’‘f does not transform nicely under the residual time diffeomorphisms, because the background 
does not. 
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(i O j, {ij) O (A;/), and so on). Since the factors of vr*-^ they are contracted with have 
precisely the same permntation symmetries, we can simply write 




= Fo ■ [l + |e7r** + a4^{5ij5ki + 

+ «6 + P2S^,SkmSln + + ■■■], 


where the a’s and /d’s are generic dimensionless coefficients, with a weak time-dependence 
that can be neglected to lowest order in slow-roll. 

If for the moment we ignore the metric perturbations 5N and iV®, then vr®-^ is simply 

TT®-^ = (cl®7r-^ -|- 5%®) -|- dkTT^dkTT^ — , (4.8) 


and isolating the different orders in vf in the action above is immediate. At the quadratic 
level we get 


C 2 


where 


Fq ■ 1^ — tt ^ -|- 2q;4/3i) idiTTj) -\- 2ci4(2 -f j3i) [diF)‘^ j 

2 ,2 2 N 


-|ea^Fo 


TT 


C 


2 

T 


C 


2 

L 


^ 6a4/3i 

e 

^ ^ 1204(1 -|- /3i) 
e 


(4.9) 


(4.10) 

(4.11) 


are the transverse and longitudinal phonon speeds. For these to be between 0 and 1, we need 
both 04 and 04 / 3 i to be small, 

04 , 04 /?! = 0{e) , (4-12) 

in analogy with the Fy + Fz = 0{e) ■ F requirement of the original solid inflation case [1]. 
In fact, in the Appendix we prove that the two propagation speed are related by the same 
constraint as in solid inflation, eq. ( 2 . 12 ), so that—as anticipated—up to quadratic order in 
perturbations our model is indistinguishable from solid inflation. 

Eq. (4.12) implies that, to lowest order in slow-roll, only the second line in (4.7) contributes 
to the cubic Lagrangian: 


£3 ~ aeFo ■ 


(8 - /^s) {diTT^Y + 4/^2 {diF) + (4/?2 + S/^s) diFdjir^dkTr^ 


G/SsdjFdjTT^dkF + 8/?4 


(4.13) 


But we are not done yet. Ref. [9] argued that the approximate internal scale invariance (2.7) 
manifests itself on the structure of the phonon self-interactions in the following way; the 
cubic action expanded about a phonon background vfo cannot correct the quadratic action 
for the fluctuations if the background is isotropic, clj 7 r° oc 6ij. Applying this requirement to 
our cubic action yields two constraints on the /?’s, 

72 + 28/?2 + 12/?3 + 48(7 + 2)/?4 = 

12/?2 + 12/?3 + 24(7 + 2)/?4 = 


11 


0 

0 , 


( 4 . 14 ) 

( 4 . 15 ) 










which allow us to eliminate (32 and 13i, 

[3a 


3/^3 - 18 
9 - 2/^3 


(4.16) 

(4.17) 


7 + 2 


We are thus left with only two free coefficients, ag and /dg, which from now on we will 
simply call a and (3. In conclusion, the cubic Lagrangian for icosahedral inflation reads 


£3 = a+o ■ {S — (3) [diTi^Y + {12(3 — 72) diTi^ [djTi’^Y + {15(3 — 72) di'Ti^djTi’^dkTi^ 



(4.18) 


Recall that in the original solid inflation model there was only one free coefficient at this 
order, Fy, appearing as an overall factor in front the cubic Lagrangian (2.19)—all the relative 
coefficients of the different terms were completely fixed. If we set our anisotropic structure to 
zero by setting (3 = 9/2, we recover precisely those ratios, and we get 


2 Fy 


W = 9 / 2 ). 


( 4 . 19 ) 


a = 


On the other hand, we will show in the next section that the choice (3 = 8 characterizes the 
completely anisotropic case, in the sense that the resulting three-point function has exactly 
zero overlap with all those that one could get from isotropic models. The cubic Lagrangian 
in this case is 


£3 = 8aFo ■ 3c?i7r* + 6 ding'djTi’^dkTi^ 



5 The size and shape of non-gaussianities 


Like in the original case of solid inflation, and for the same reasons spelled out there [1], the 


leading trilinear interactions we need to consider to compute the scalar three-point function 
are the phonon self-interactions we wrote down above. That is, we can neglect interactions 
involving the metric perturbations. The computation of the three-point function parallels 
that in [ 1 ], with obvious modihcations due the new tensor structures we have in the cubic 
Lagrangian. Neglecting the weak time-dependence of the scalar modes outside the horizon, 
the result is 


(C1C2C3) — { 2 Tr)^ 6 ^{ki + ^2 + ^3) X 

9 a Q{ki,k 2 ,kz)U{ki,k 2 ,h) 


(5.1) 


32 M4 ■ e^cf ' kfkiki 
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where 


Q{ki,k2,h) = (8 -/3)/cifc 2 A ;3 + 6/3 


[ki ■ ^ 2 ) {k2 ■ ks) (ks ■ ki) 


{h ■ hY 

k2k‘i 


-\- h‘ 


kihh 
{h ■ kiY 


+ k 


(ki ■ hY 


kski kik2 

+2)k\k\kikikiki + ( 7 + i){k\k\k^kiki+^ki+^ 


+ (9/3 - 48) (k 

-16/3 + 72 1 

7 + 2 kik2k:i 

P P -I- P P P 4- Ah^ P 

-I-4P P P P -I- 4P P+1 P t*+l P1 -L It® t® P P P“1 P“1 
/i/ 2 ^ 2 ^ 3^3 a^2^2 ^ 3 ^ 3 / \^l^l^ 2 ^ 2^3 ^3 

_l_Li ui ^®-lp-lp p _|_ ^®-l^®-l^® P P P _|_ At® t® t® P“1P P“1 

A/A^A/2 ^2 ^3^3 ^2^2^3^3 A/2^2 ^3^3 

+Ak\k\-^kikikiki-^ + Ak\k\-^k\k\-^kiki) 


and 


U{ki,k2,k3) = 


3{kl + kl + kl) +20klklk. 


/Ci/C2^3 (^1 + ^2 + ^ 3 ) 

+18(fc^/c2fc3 + kik^ks + A;iA:2fc3) + 12(A:^fc2 + ^2^3 + ^3^1) 
9 (A:^fc 2 + Sperms) + 12 (A:^A :2 + Sperms) 

+18(fc^/c2A;3 + Sperms) I . 


(5.2) 


(5.3) 


The overall delta function leaves us with only two independent momenta, say k 2 and k^. 
Usually, because of isotropy, the absolute orientation of these two vectors does not matter, 
and one needs only three independent quantities to characterize the kinematical configuration: 
the magnitudes k 2 and /C 3 , and the relative angle 6 . For us, because of our anisotropies, 
the absolute orientation matters, and so we have to keep all the six components of k 2 and 
k^. This complicates the analysis considerably. In particular, we cannot use the standard 
techniques of [ 6 ]. 

A convenient parametrization of ^2 and k^ is the following one. Dehne 62 and 02 as the 
standard polar and azimuthal angles of /c 2 ; but define 6*3 and 03 as the polar and azimuthal 
angles of k^ with respect to a primed coordinate system in which the z' axis is along ^ 2 , 
and the x' axis lies in the plane dehned by the 4 and z' axes (see hg. 2). The cartesian 
components of k 2 and k^ thus are 


^2 = ^2 (sin 6*2 cos 02 , sin 6*2 sin 02 , cos 6 * 2 ) (5.4) 

^3 = ^3 (sin 6*2 cos 02 cos 6*3 + cos 6*2 cos 02 sin 6*3 cos 03 — sin 02 sin 6*3 sin 03 , 
sin 62 sin 02 cos 6*3 + cos 6*2 sin 02 sin 6*3 cos 03 + cos 02 sin 6*3 sin 03 , 
cos 62 cos 6^3 — sin 62 sin 6*3 cos 03) . (5.5) 


The advantage of this parametrization is that 62 , is the relative angle between ^2 and k^, and 
so any dependence on 6*3 is perfectly consistent with isotropy. Anisotropies show up as a 
non-trivial dependence on 02 , 03 , and 6 * 2 . 
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Figure 2: The coordinate system defined in the text. The z, z' = k 2 , and x' axes all lie in the same 
plane. 


Following the standard conventions for correlation functions of the Newtonian potential 


$ = 

ic 

(5.6) 

($(fc;)<i>(fc;)> = 

(27r)^(5^(fci ^ 2 )^ 

(5.7) 

($(fc;)<i>(/r2)$(/r3)> = 

(27r)^(5^(fci + k2 + k2)f{ki, 4, h) , 

(5.8) 


we dehne /nl nsing eqnilateral conhgnrations: 

^ ^ ^ 6A^ 

= hL-^ . (5.9) 

However, the equilateral-triangle condition only fixes the relative angle 63 , and so in our case 
the resulting /nl depends non-trivially on the other angles, 02 , 03, and 6 * 2 . To get a readable 
expression, we average /nl over 02 and 03 , 

/nl(6 *2) = y d02(i03 /nl(6 '2, 02, 03) • (5.10) 

The remaining dependence on 62 will still be a measure of anisotropy. For our three-point 
fnnction, we get: 


A(j, 


f{ki,k2,h) 


/nl(6'2) 


9 1 

Too ^ 

15 a 2 Q{ki,k 2 ,kfijU{kx,k 2 ,kfis 

o ,„2 ■ ■ A’3A’3 l3 

+w(2/^-9)^6(cos02)] , 


(5.11) 

(5.12) 

(5.13) 
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where Pg is the sixth order Legendre polynomial. 

The typical size of /nl is the same as in the standard solid inflation case, parametrically 
as big as l/ec| if one assumes a ~ 1 (analogous to Py ~ F for solid inflation). But clearly 
the most interesting feature here is the angular dependence of /nl: the appearance of Pg in 
/nl(^ 2 ) is the first indication that the case with /3 = 8 is a very special one, with a completely 
anisotropic /nl: if we average /nl over cos 02 , which is equivalent to averaging the full /nl 
over all angular variables, we get zero. 

We can go further and, following [6], consider the overlap between our three-point function 
and other ‘shapes’. This is defined as 

oos(/./')-^=T^ (5.14) 


where 

f ■ f = ^f{ki,k2,h)f'{ki,k2,k3)/{al^al^alj . (5.15) 

ki 


The sum runs over all triangles in momentum space, and is in fact an integral since the 
momenta are continuous variables. 

If for / we take our shape (ignoring overall constant factors, which do not contribute to 
the overlap (5.14)), 


/(fci, ^2, ^ 3 ) 


Q(fci, k2, h)U{ki, k2, ks) 

1^0 T^o T^o 
^ 1 ^ 2 ^3 


(5.16) 


and for /' that coming from a general isotropic model. 


f f{ki, k 2 , ks) , 


(5.17) 


we hnd exactly vanishing overlap if /3 = 8. Again, the reason is manifest if, when computing 
the angular integrals for the overlap (5.14), we perform the integrals over 02 and 03 hrst: 


J dcos 02 d 02 dcos 03 003 f f' 


a S27i'^U{ki,k2,ks) 

dcos 63 — - di d 7 d —^^ 2^3 f (ki, ^ 2 , k^) 


(5.18) 


7 

I A/2 

X / 0COS 02((0 — 8)Gi(A:2,/::35 02, 03 ) + (20 — 9)(j'2(^2, ^3; 02, 03 )) , 


where 


Gi{k 2 , k^, 02 , 03 ) 
^ 2 (^ 2 , ^3, 02 , 03 ) 


— [2(^2 + ^3) P2(cOS 03) -|- ^2^3 cos 03(1 -|- 3 COS^ 03)] (5.19) 

[^2-P2(cOS 03) 2fc2/C3 P3 (cOS 03) fcg P4 (cOS 03)] ■ Pg(cOS 02)(5.2O) 


and the P„’s are the Legendre polynomials. All quantities that only depend on the magnitudes 
ki, k 2 , ^3 factor out of the 02 integral, because they cannot depend on the orientation of the 
triangle defined by the momenta. 

We clearly see the two limiting cases now. For 0 = 9/2, only the Gi contribution survives, 
it has no 02 dependence, and we recover the results of the isotropic solid inflation case. On 
the other hand, for 0 = 8, Gi is gone and G 2 , being proportional to Pq{cos 92 ), averages to 
zero when we integrate over 02 . 
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6 Anisotropic tensor spectrum? 

The existence of an anisotropic six-index invariant tensor snggests that anisotropies can 
also show np in the tensor modes’ two-point function, because of the possible quadratic 
Lagrangian term 

( 6 . 1 ) 

However, it is easy to convince oneself that such a term cannot arise from expanding the 
lowest-derivative action we have been working with so far, 

S = J d^x^[^M^R + F{B^^)] , (6.2) 

simply because all possible anisotropies are in the structure of F, but its argument = 
dcj)'^ does not involve derivatives of the metric. 

On the other hand, in the presence of higher derivative terms, one will generically get 
such a term. Consider for instance the invariant 

(9 • • • , (6.3) 

where {g ■ ■ ■ g) stands schematically for any twelve-index tensor built out of the metric. Setting 
the 0^’s to their background values , and expanding in powers of the tensor modes 7 , the 
covariant derivatives have the schematic form 

VV0^ ^ H + d-f; (6.4) 

and so, upon taking all the contractions in (6.3), one does expect to hnd the term (6.1) at 
quadratic order. Similar considerations apply to higher-derivative terms that involve higher 
powers of curvature tensors, for instance a trilinear term schematically of the form 

d^ct>^d,ct>^dpct>^d„ct>^f , (6.5) 

with suitable contractions with our anisotropic invariant tensor T/l'), '/. (We need at least 
three Riemann tensors, because our Taniso is totally symmetric, while has antisymmetry 

properties as well.) 

However, if we want the anisotropic quadratic terms that we get from these higher 
derivative corrections to compete with the purely isotropic ones we get from the Einstein- 
Hilbert action, we need to give the higher derivative corrections a large coefficient, of order 
Mp/H^ in the examples above. This makes the smallness of F, 

F ~ , (6.6) 

potentially unstable against quantum corrections. For instance, we expect graviton loops 
involving 0’s on the external legs and the coupling (6.5) in the vertices, to drastically correct 
F[B^^). A quick order-of-magnitude estimate of this two-loop diagram 


gives a correction to the effective action 

AF ~ ^£21/2 ^ 

H 

assuming the phonon speeds are relativistic, cl,ct ~ 1, and we cutoff the loop integrals at 
the solid’s strong coupling scale, Agtrong ~ [1], For AF to be at most of order F, we 

need a small enough e: 

£ < ■ ( 6 . 8 ) 

This goes in the opposite direction to the bound on e that guarantees that perturbations are 
weakly coupled at freeze out [1], 

(6.9) 

but, for small if/Mp, it is perfectly compatible with it. 

We will not attempt a systematic analysis of higher derivative corrections here. However, 
the above estimates suggest that it may be consistent to expect higher derivative corrections 
to be large enough to yield order-one anisotropies in the tensor spectrum, but small enough 
to preserve the technical naturalness of our effective field theory. We will investigate the 
issue further in a forthcoming publication, where we will also give the detailed form of the 
resulting anisotropic tensor spectrum. 

7 Concluding remarks 

We have shown that, among the generalizations of solid inflation with discrete rotational 
symmetries, that with icosahedral symmetry is the only one that is naturally compatible 
with the observed isotropy of the background and of the scalar spectrum. 

The associated scalar three-point function is in general highly anisotropic, and this 
suppresses its overlap with all the standard templates used in CMB data analyses. For 
a specific choice of the Lagrangian coefficients = 8, in our notation), it is completely 
anisotropic, in the sense that such an overlap vanishes exactly. This leaves open the possibility 
that a large non-gaussian signal is hiding in the data, waiting to be unveiled by a dedicated 
anisotropic analysis. 

It is worth pointing out that our anisotropies are not of the same nature as those discussed 
in [10, 11]: there, for any given realization one expects small anisotropies in the scalar 
spectrum; but, statistically speaking, these average to zero. On the other hand, in our case it 
is the statistical correlation functions themselves that are intrinsically anisotropic, potentially 
maximally so. 

Similar considerations apply to the tensor spectrum as well: in the presence of sizable 
higher-derivative corrections, it can be highly anisotropic, which makes the standard detection 
strategies inefficient, and calls for a dedicated analysis. 

For the scalar modes, we see no reason for the /3 = 8 case to be preferred over others; 
for instance we see no symmetry protecting it against quantum corrections. However, it is 
a simple, consistent limit of our theory, and we find it interesting that such a completely 
anisotropic limit exists at all. Is it an accidental feature of our truncation of the theory at 
the cubic/three-point function level, or does it survive at higher orders as well? 
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Appendix 

A Background stress-energy tensor 

To relate the background stress-energy tensor to F and its derivatives, it is useful to organize 
F’s dependence on in terms of the variables 

nij 

X = IB], = (B" = !f‘-d,4>‘d„4>-’)- (A.l) 

X depends on the overall normalization of whereas 6^*^ does not. As a result, the 
approximate internal scale invariance (2.7) translates into a weak X-dependence of F. 
Taking the variation with respect to the metric for the solid action 

S-soiid = j d^x^F{X, , (A.2) 

we hnd the stress-energy tensor 

= g^uF - 2Fxd^(t)^du(t)^ - ^Fij[d^(t)^dyct)-^ - 8^(1)^) , (A.3) 

where the subscript X and IJ stand for partial derivatives w.r.t. X and b^^. 

When we evaluate on the background conhguration, we can use the fact that Fjj must 
be icosahedral invariant. As we saw, for a two-index tensor this implies that it is proportional 
to Sjj. The terms in parentheses in (A.3) thus cancel against each other, and we are left with 
the same background stress-energy tensor as in 5'0(3)-invariant solid inflation. In particular: 

p = -F, p = F- ^Fx 
FxX 


B Phonon propagation speeds 

To hnd the phonon propagation speeds, we should expand the solid action (A.2) to quadratic 
order in the phonon held tt. Let’s use the same X and variables of last section; the 
expansion of the action then reads 

C = Fx6X + F,j 5F^ + iFxxiSXf + Fx,^J 5X 6F^ + \F,,^ki SF^Sb^^ + ... (B.l) 


(A.4) 

(A.5) 
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When we specialize all the derivatives of F to the background, by icosahedral symmetry they 
must take the form 


^iji ^X,ij ^ij ) Pij,kl fl “ 1 “ ^2 F ^il^jk^ , 


(B.2) 


with generic, time-dependent coefficients. This kills some of the terms in (B.l) because, by 
definition (eq. (A.l)), the fluctuation of W is traceless. We are left with 


C ~ Fx SX + \Fxx(6XY + h (»■')" 


(B.3) 


We thus need 5X up to quadratic order in the phonon field, and up to linear order. 
These are 

hX = TT** , 




TT*-^ = (9V-^ -I- -1- dnTT^d^TT^ 


where is the fluctuation of 

+ TT*^' , 

At quadratic order in n we get 

C 2 = -Fx[7r^ - CT^diTijf - {cl - 4)(V ■ vf)^ ] 

with 

-L 


9 Fxx 8/2 9 2/2 

c?=l + 2^^ + —— 4 = 1 + - — 

Fx 27Fx’ ^ 9Fx 


This is identical to the S'0(3)-invariant solid inflation’s result, upon identifying 


/2 


here 


{Fy + Fz) 


there 


(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B,8) 


As a consequence, icosahedral inflation still obeys the relation (2.12), and, more in general, is 
indistinguishable from solid inflation at the quadratic level. 
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